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Abstract 

Axisymmetric equilibria with incompressible flows of arbitrary direction are 
studied in the framework of magnetohydrodynamics under a variety of physi- 
cally relevant side conditions. To this end a set of pertinent non-linear ODEs 
are transformed to quasilinear ones and the respective initial value problem is 
solved numerically with appropriately determined initial values near the mag- 
netic axis. Several equilibria are then constructed surface by surface. The 
non field aligned flow results in novel configurations with a single magnetic 
axis, toroidal shell configurations in which the plasma is confined within a 
couple of magnetic surfaces and double shell-like configurations. In addition, 
the flow affects the elongation and triangularity of the magnetic surfaces. 
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1. Introduction 



In a previous paper [T] the first two authors derived a generalized Grad- 
Shafranov equation governing the magnetohydrodynamic (MHD) equilibrium 
states of an axisymmetric plasma with incompressible flows [Eq. ([1]) in sec- 
tion 2]. By assignment of the free functions contained in (pQ), known solu- 
tions of the Grad-Shafranov equation, i.e the Solovev solution [2} [3] and the 
Hernegger-Maschke solution [U [5] , were extended in Refs. [6] and [7] , respec- 
tively. For both extended solutions the flow can change the magnetic field 
topology, thus resulting in a variety of new configurations of astrophysical 
and laboratory concern. 

Instead of assigning the free surface functions of (pQ), it may be of physical 
or mathematical importance to introduce side conditions, e.g. isodynamicity: 
B 2 = B 2 (ip), where if) and B are the poloidal magnetic flux function and the 
magnetic field modulus, respectively. In the quasi-static case, viz. when the 
flow is neglected in the momentum equation but it is kept in Ohm's law, 
it was proved in Ref. [8] that there is a unique configuration of this kind 
with circular magnetic surface cross-sections near the magnetic axis. This 
equilibrium was fully constructed in Ref. [9]. The same configuration persists 
in the case of flows aligned to the magnetic field [1] . In the case of non field 
aligned flows satisfying the side condition P + B 2 /2 = f(ip), where P is the 
thermal pressure and / is an arbitrary smooth function of ip, the magnetic 
surfaces near axis become elliptical with elongation perpendicular to the axis 
of symmetry [TU] . 

The aim of the present study is to construct axisymmetric steady states 
with incompressible flows of generic direction under a variety of side condi- 
tions and to examine the impact of the flow on the equilibrium characteristics 
and particularly in connection with the magnetic topology. A preliminary 
investigation was conducted in Ref. [UJ. Here, the construction is carried 
out numerically on the basis of a procedure suggested in Refs. [8], [9J and pQ. 
The main conclusion is that the flow results in a variety of novel equilibria 
and opens up the possibility of changing the magnetic field topology. 

The side conditioned equilibrium equations are briefly reviewed in section 
2 along with the solving procedure of Refs. [S] and pQ. In addition, the orig- 
inal ODEs of concern are mapped to quasi-linear ones by a transformation 
which for the quasi-static case was employed in Ref. |9j. After establishing 
initial values of the unknown functions near the magnetic axis the problem 
becomes well posed and is solved numerically. The various kinds of config- 
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urations associated with numerical solutions are presented in section 3 and 
are compared with the existing ones in the literature. Section 4 summarizes 
the study and the conclusions. 

2. Side conditioned equilibria 

The MHD equilibrium states of an axisymmetric magnetized plasma with in- 
compressible flows are determined by the generalized Grad-Shafranov equa- 
tion pQ, 



s + 2 VI — M 2 



along with the Bernoulli relation for the pressure 

P = P S (4>)- P 



2 R 2^2 



1 - M 2 



= 0, 
(1) 

(2) 



Here, (z, R, <ft) are cylindrical coordinates with z corresponding to the axis 
of symmetry; the function ip(R, z) labels the magnetic surfaces; M{ip) is the 
Mach function of the poloidal velocity with respect to the poloidal-magnetic- 
field Alfven velocity; p(V') and are the density and the electrostatic 

potential; X(ijj) relates to the toroidal magnetic field; for vanishing flow the 
surface function P s (ip) coincides with the pressure; v is the velocity modulus 
which can be expressed in terms of surface functions and R; A* = R 2 V ■ 
(V/-R 2 ); and the prime denotes a derivative with respect to ip. Derivation 
of ([1]) and ([2]) is provided in Ref. pQ. The surface quantities M(ip), 
X(ip), p(ip) and Ps{tp) are free functions for each choice of which ([T|) is fully 
determined and can be solved whence the boundary condition for ip is given. 
By using the transformation 
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JTJ) and (T5]), respectively, reduce to 
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Note that no quadratic term as |Vw| 2 appears anymore in (111) . The forms of 



and §5§ indicate one to introduce the new surface quantities 



N(u) 



X 



VI - M 2 ' 



L(u) = VP 



du ' 



(6) 



which are helpful in reducing the number of the explicit free functions by 



one. 



Instead of specifying the free functions to determine (pQ), one can introduce 
side conditions as those quoted in table [TJ They can be expressed in terms of 
the thermal pressure P, the magnetic pressure B 2 /2, the flow energy density 
pv 2 /2 or combinations of them and consist in that these quantities remain 
uniform on magnetic surfaces. Such conditions lead, in general, to an 
additional relation between (Vu) 2 ,u and R as already accomplished in Ref. 
P]. Specifically, (jlj) and (jSJ) can be put in the respective forms 



|Vm| 2 = 2[i(u) + R 2 j{u) + R 4 k(u)} 
A*u = -f(u) - R 2 g{u) - R 4 h{u) 



where 
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(7) 
(8) 

(9) 
(10) 



(12) 



The other coefficients j(u) and k(u) being side condition dependent, are given 
in table [H 

Equations © and (|SJ) can be solved simultaneously by the method sug- 
gested by Palumbo [5] according to which the variable z is treated as function 
of u and R. The method leads to compatibility conditions in the form of five 
ordinary differential equations (ODEs) [Eqs. (44)- (48) of Ref. jTj] along with 
an integral relation for the function z(u,R) [Eq. (49) of Ref. p]]. Thus, to 
determine the equilibrium a) the ODEs should be solved on each magnetic 
surface and then b) the intersection points of that surface with the poloidal 
cross-section should be obtained by the integral relation. There is a difficulty, 
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however, stemming from the non linearity of the ODEs. Here, to overcome 
this difficulty, the function u is mapped to a new function w through the 
transformation: 

''" -( 9+ *W„), (13) 



dw \ du r 

where J~{w) is an arbitrary smooth function. 

Under this transformation the compatibility condition leads to the fol- 
lowing set of ODEs 

4Q + 2QY' - Y(Y + &') = 0, (14) 

82 + F(4 - Q') + 2QY' + 6' = 0, (15) 

- 3 - 2&H' + Q' + 2EY' + H(QY + &' - YE') = 0, (16) 

- 2QH' + H(Q' - 8) + E' = 0, (17) 

- 5H 2 - 2EH' + HE' = 0, (18) 

where the functions Q(w), Q(w), Y(w), E(w) and H(w) can be expressed in 
terms of the initial physical surface quantities. Note that ffT4"|) - ffTBl are quasi- 
linear, viz. the derivatives appear linearly, and therefore Picard' s theorem 
guarantees existence and uniqueness of the respective initial value problem. 
Because of indefmiteness of some of (|T4 l) -f |T8l) on the magnetic axis when 
they are put in solved forms, initial values of the unknown functions near 
the magnetic axis can be obtained on the basis of Mercier expansions around 
axis along with a l'Hospital-like procedure. Details will be given elsewhere. 
Also, the integral relation for z(w,x) is given by 

p ±j[Hx 2 + x-Y] 



dz 
dx 



V [& + Qx + x 2 E - \x (Hx 2 + x - Y) 2 ] 



21 1/2 ' 



(19) 



where x = R 2 , p = dw/dx, q = dw/dz, r = d 2 w/dx 2 and q = d 2 w/dz 2 . 
Once initial values are established the problem is well posed and can be solved 
numerically. Accordingly, we have developed a programme in Mathematica 
5.1. For up-down symmetric configurations to be considered here there are 
three free parameters (Rq, E and Ho) associated with the radial distance of 
the magnetic axis and the functions E and H thereon. It is also noted that 
the problem can be solved by the transformation, alternative to ( TL3l) 



du 1 /, dk\ _ . . . 

--(h+-)=F(w). (20) 



dw 2 \ du J 
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3. Magnetic configurations and impact of flow 

For the solutions to be presented here we have chosen, without loss of 
generality, R = 1 and w = 0, where w is the value of w on axis. Depending 
on the sign of S , there are two kinds of configurations: 

1. Toroidal configurations for So < having a single magnetic axis [lo- 
cated on [z = 0, R = R ) where w = w = 0]. These configurations 
have magnetic field topology similar to those of Refs. [9] and [ID] , w- 
contours of such a solution are shown in Fig. HJ In general for S < 
the magnetic surfaces are more elongated horizontally (parallel to the 
mid-plane z = 0) up to the magnetic compared with the quasi- 
static ones corresponding to S = H = 0. 

2. Toroidal shells for Ho > in which the plasma is contained within 
two toroidal surfaces, w-contours for a solution of this kind are shown 
in Fig. [2j Although mathematically the solution has an extremum 
on w = 0, the physically acceptable part is restricted to non positive 
values of 0(u>), thus resulting in a toroidal vertically elongated shell. 

The various kinds of configurations can also be classified in terms of the 
parameters S and H as follows. 

1. S = and H ^ 0: 

Let us first note that S(w) = implies parallel flows (or the quasi- 
static equilibrium) because then it follows from ffTS]) that H(w) = 0. 
This should not be confused with the case of H = (on axis) in this 
paragraph which involves non parallel flows when H ^ 0. As in the 
case of parallel flows, however, the magnetic surfaces near axis have 
circular cross sections. For H > the surfaces far from axis are less 
parallel elongated than those of the quasi-static equilibrium. For H < 
the triangularity can change drastically. As an example, a configuration 
with inverse triangularity is shown in Fig. [31 

2. S ^ and H = 0: 

Equilibria for So < and Hq = have been constructed in Ref. [10] by 
a different method. Unlike in Ref. [10J, however, no restriction on the 
elongation of the magnetic surfaces (parallel to the mid-plane z = 0) 
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was found here. As a matter of fact such a configuration with very 
elongated surfaces is presented in Fig. @] (in purple). Thus, the limi- 
tation on the elongation reported in [10] may be due to the particular 
method of solution in that paper. For So > the equilibrium becomes 
a toroidal shell with magnetic surfaces elongated perpendicular to the 
mid-plane z — 0. 

3. S ^ and H ^ 0: 

In this generic case it is particularly interesting to examine whether 
there are configurations with two magnetic axes. This requires two 
roots in the numerator of (fT§j) and four roots in the denominator ap- 
propriately located with respect to the roots of the numerator. (Note 
that for a quasi-static equilibrium (H = E = 0) only configurations 
with a single magnetic axis are possible.) For this reason we first ex- 
amined this requirement by applying the Sturm theorem and Descartes 
rule [12]. The former determines the exact number of real roots of a 
polynomial with real coefficients; the latter determines the maximum 
number of positive roots of such a polynomial. It turns out that the 
requirement is compatible with the Sturm theorem and Descartes rule. 
Then, by inspection we found that the requirement can be fulfilled for 
S > and H < 0. An equilibrium of this kind shown in Fig. con- 
sists of a toroidal shell reaching the axis of symmetry, similar to those 
reported in paragraph 3.2, and a second thin shell-like configuration lo- 
cated farther from the axis of symmetry. The distance between the two 
configurations decreases as \E /H \ takes larger values. As can be seen 
in Fig. 0, however, the magnetic surfaces of the shell-like configuration 
do not close. Closeness does not improve either by varying So and Hq 
or by using a different numerical method in FORTRAN. Therefore, the 
existence of double toroidal shell configurations with closed magnetic 
surfaces remains an open question. 

For the other three combinations of signs of S and H one can obtain 
configurations similar to those presented in paragraph 3.2 As in the 
case of H = 0, for S < there is no limitation on the elongation of 
the magnetic surfaces parallel to the mid-plane z = as |So| increases. 
Also, it is worth to mention the strong change in the triangularity of 
the configuration for Hq < as it is illustrated in Fig. [6j 

4. Summary and Conclusions 
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We have studied axisymmetric equilibria with incompressible flows under 
side conditions of physical relevance by a procedure introduced in Refs. [8], 
[9] and [JJ . This procedure reduces the problem to a set of ODEs for certain 
surface functions and an integral relation determining the points of the cross 
section of a magnetic surface with the poloidal plane. Because of the nonlin- 
earity of the original ODEs, we have employed transformation f[T3"j) mapping 
the original ODEs to quasilinear ones [Eqs. (|14l) - (|18p ] and containing equal 
number of unknown surface functions; thus, existence and uniqueness of the 
respective initial value problem is guaranteed. After determining appropri- 
ately initial values near axis, because of indefiniteness of the ODEs thereon, 
the problem has been solved numerically surface by surface with two free 
parameters (H and H ) associated with the non- field aligned flow. 

The flow results in the following novel kinds of up-down symmetric equi- 
libria: 

1. Configurations with a single magnetic axis for S < with magnetic 
surfaces near axis elongated parallel to the mid-plane z — 0. For So = 
the magnetic surfaces near axis become circular as the quasi-static 
ones corresponding to H = H = 0. The special case of equilibria 
constructed by a different method in Ref. [TU] are recovered for H = 0. 
Unlike in Ref. [TU], however, no restriction on the elongation has been 
found in the present study. 

2. Toroidal shells for S > and H > in which the plasma is confined 
in the interior of two nested magnetic surfaces. The magnetic surfaces 
are elongated perpendicular to the mid-plane z = compared with the 
quasi-static ones. 

3. Equilibria consisting of a toroidal shell reaching the axis of symmetry 
and a second shell-like configuration for S > and H < 0. Thus, the 
flow opens up the possibility of changing the magnetic field topology. 

Also, the shape of the magnetic surfaces far from axis is affected by the value 
and sign of H ; specifically: a) they become less elongated parallel to the 
mid-plane z = as Hq takes larger positive values and b) the triangularity 
of those surfaces is affected drastically for negative values of H . 

It is emphasized that the above reported conclusions hold irrespective 
of the particular condition of table [1] except for P + B 2 /2 being uniform on 
surfaces which corresponds to Hq = 0. The properties of particular equilibria 
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in connection with profiles of the (original) physical quantities, i.e. pressure, 
magnetic field, velocity etc, deserves further investigation. Also, in view of 
the tough and in general unsolved stability problem of steady states with flow, 
the stability of the equilibria constructed here remains an open question. 
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Table 1: The coefficients j(u) and k(u) for several side conditions involving 
the quantities P, B 2 and pv 2 . Note that the condition P + pv 2 /2 = f(u) 
not included in the Table can be satisfied only for parallel flows [cf$ jdu = 0] 
because of the explicit independence of the last term in ([5]) . 
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Figure 1: Two sets of w-contours the one (in yellow) for S = —0.5 and 
H Q = 0.5 and the other (in purple) for H = and H = 0.5. Both equilibria 
have a magnetic axis on which w — 0. 
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Figure 2: A set of w-contours associated with a toroidal shell for S = 0.1 and 
H Q = 0.5 (in yellow). The second set of u>-contours (in purple) corresponds 
to an equilibrium with a magnetic axis (S = and H = 0.5). 

z 
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Figure 3: A configuration with inverse triangularity (in yellow) for S = 
and H Q = —0.4. The u>-curves in purple correspond to the quasi-static 
configuration. 
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Figure 4: Two configurations with a magnetic axis and very elongated mag- 
netic surfaces parallel to the mid-plane z = 0, the one for S = —10 and 
H = (in purple) and the other for S = — 10 and H = 1 (in yellow). The 
respective values of the eccentricity are 0.988 and 0.9. 
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Figure 5: An equilibrium for S = 0.05 and H = —0.1 consisting of a 
toroidal shell and a smaller toroidal shell-like configuration. 
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Figure 6: Impact of H < on the triangularity of an equilibrium with a 
magnetic axis (S = —0.1): a) a quasi-static like configuration (H = —0.09); 
b) an eye-like configuration (H = —0.3); and c) a configuration with inverse 
triangularity (H = —0.4). 
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